A generalized Coulomb logarithm for strongly coupled plasmas 
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Transport coefficients for Coulomb collision processes, e.g., friction, energy exchange, and resistivity, are cal- 
culated for Debye-screened (Yukawa) plasmas including correlation effects within the binary collision approxi- 
mation. A generalized Coulomb logarithm, 5(A), is required, which asymptotes to In A in the weakly coupled 
limit, but is significantly modified when A < 10 2 due to large angle collisions. In the regime 10~ 2 < A < 10 2 , 
S can be substantially larger for attractive collisions than repulsive. 
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Plasmas in several modern research areas, including 
inertial confinement fusion (ICF)P antimatter plasmas, 2 
ultracold plasmas,^ and dusty plasmas, 4 can enter 
regimes in which the central tenet of conventional plasma 
kinetic theory is violated. Conventional theories^"^ are 
based upon a weak coupling approximation, often re- 
ferred to as "the plasma approximation," which assumes 
In A » 1, where A ~ 12t™A^ ~ T 3 / 2 /(Z iy /n) is the 
plasma parameter, and Ap = E s ^e> 2 ) ^ 2 ^ s ^ ne ^e- 
bye length (A 2 ^ = T s /4Trq^n s ). Most plasmas in space 
and the laboratory are weakly coupled, In A > 10, but 
moderate, 2 < In A < 10, and strong, In A < 2, cou- 
pling regimes can be reached if the plasma is extremely 
dense (ICF), cold (antimatter and ultracold), or highly 
charged (dusty). Understanding how Coulomb collisions 
affect transport properties, such as friction, temperature 
equilibration, and electrical conductivity, is essential in 
order to accurately predict the behavior of these plasmas. 

A generalized method of calculating transport coef- 
ficients for fluid equations, which relaxes the conven- 
tional weak coupling approximation, is provided. The 
weak coupling approximation is utilized twice in con- 
ventional theories, and both instances are avoided here. 
The first instance is that conventional theories are based 
on a small scattering angle expansion, or a hierarchy 
expansion (BBGKY), which orders terms according to 
In A 3> 1. The second instance is that a logarithmi- 
cally divergent integral arises as a result of using the 
bare Coulomb potential for individual particles. Debye 
screening is introduced, and the integral kept finite, by 
imposing an ad-hoc cutoff for the maximum impact pa- 
rameter at the Debye length. This cutoff is only justi- 
fied when In A 3> 1. Here, we avoid the small scattering 
angle expansion by exploiting symmetries that arise in 
velocity- space moments of the full Boltzmann collision 
integral. The theory self-consistently accounts for Debye 
screening, and thus avoids divergent integrals, by using 
the screened Coulomb potential. 

The screened Coulomb potential is common, but not 
universal, in the plasmas mentioned above. For example, 
collective interactions can cause more complex behavior 
than Debye screening, such as wakes surrounding fast- 
flowing particles. 81 These effects are especially important 
for fast particle stopping power in ICF. However, it is a 



valid approximation for many processes of interest, re- 
gardless of correlation strength. 

We find that a generalized Coulomb logarithm, S, is 
required. This asymptotes to 2 ~ In A in the weakly 
coupled limit, but is modified significantly by correlation 
effects when A < 10 2 . In the regime, 10~ 2 < A < 10 2 , 3 
can be substantially larger for attractive collisions than 
repulsive. The present work focuses on species with flow- 
ing Maxwellian distributions. However, the underlying 
approach does not depend on the functional form of the 
distributions and can be applied more generally. One 
such example is Spitzer's approach^ which considers de- 
viations from Maxwellian as higher-order perturbations. 

The theory is based on the Boltzmann collision op- 
erator, which requires that the plasma be sufficiently di- 
lute to justify the binary collision approximation. Effects 
of three, and more, particle collisions, as well as phase 
transitions, become important at sufficiently strong cou- 
pling. Nevertheless, this approach is valid for coupling 
parameters (A -1 ) orders of magnitude larger than con- 
ventional theory can describe. To demonstrate this, we 
compare the calculated temperature relaxation rate be- 
tween electrons and ions for ICF relevant parameters 
with recent molecular dynamics (MD) simulations!^ Our 
theory agrees with the simulation results for all values of 
A that were simulated (10 _1 < A < 10 3 ), whereas con- 
ventional theories widely disagree with them for A < 10. 
Fluid codes used to simulate strongly coupled plasmas 
require transport coefficients that can e xten d beyond the 
weakly coupled Landau-Spitzer theoryP^l Current fluid 
codes disagree with MD simulations of thermal equili- 
bration rates, hot spot formation, laser absorption, and 
neutron yield in correlated regimes; processes essential to 
accurate modeling of ICF implosions EDEU 

The Boltzmann kinetic equation describing the evo- 
lution of the distribution function of a species s is 
dt fs + v ■ V/ s + a • V v / S = J2 S > CbC/s, fs'), in which 

C B (f s Js')= [ [ dVdfitru (1) 
Jv' Jn 

is the operator describing collisions between species s 
and s'. Here f. = /,(v), f' s , = / s ,(v'), /, = / a (v), 
f s , = / s '(v'), in which (v, v') are the initial parti- 
cle velocities, and (v, v') = (v + Av,v' + Av') are 



2 



the particle velocities after a collision. The differen- 
tial scattering cross section is denoted a, the solid an- 
gle d£l — d(f>d6 sin 9, the relative velocity u = v — v', 
and u = |u|. Fluid equations of motion are derived from 
velocity-space moments of the plasma kinetic equation.^ 
Here we are interested in the effects of Coulomb col- 
lisions, which enter through the friction force density 
R s ~ s = J d 3 v m s vC(f s , f s >), energy exchange density 
Q s - S ' = jd 3 v\m s {v - V s ) 2 C(f s ,f s ,), and so on for 
higher-order moments. The reference resistivity, r/ , can 
be obtained from the electron-ion friction force density, 
R e_I 7y en e J, where J = en e (Vj — V e ) is the current 
density. 

The collision operator used to evaluate R s s and 
Q s ~ s in weakly coupled plasma theorjB^ can be ob- 
tained from Eq. ([I]) as follows. After applying a small 
scattering angle expansion (v, v') ^> (Av, Av'), integrals 
of the form 



obtains b„ 



An in the adiabatic limit, but misses b„ 



{Au} 



dfl a u Au, 



(2) 



as well as {Au. . . Au} ariseP in which Au = 
■u[sin6> cos </> x + sin # sin(/>y — 2 sin 2 (6>/2)u] is the rotation 
of the relative velocity vector after scattering through 
angle 9. Typically, the Rutherford scattering formula 
o = q 2 q 2 ,/[4m 2 sl u 4 sin i (9/2)} is used in Eq. ^ and 
Debye screening is introduced by assuming the maxi- 
mum impact parameter is the Debye length. Doing so 
yields {Au} = —Anq 2 s q 2 , hi(XY)/b min )/(m 2 s ,u 2 ), where 
^min = Qsqs' /(m SS fU 2 ). If b^ in is smaller than the de 
Broglie wavelength, 6^™ n = h/(2rn ss >u) is used. 

Using the above approximations, and keeping only 
terms of C[ln(AD/&min)], the coll ision operator can be 
written in the Landau forml^ES Taking the distribu- 
tion functions to be flowing Maxwellians, / s (v) = 
n s exp[— (v — ~V S ) 2 / v Ts]/ n3 ^ v Ts^ the friction force den- 
sity and energy exchange density are found to be 



R s 



-n s m s i/ s - s ,(V s - V s /) 



(3) 



and 



qs-s =_ 3 _=_ n4lw (T a -IV), (4) 

TO.,' 



to lowest order in In A > 1, and |V S - V s >\/v T < 1. The 
characteristic collision rate is 



16^/Ttq 2 s q 2 s , n s > „ 
3m s m ss >v^, 



(5) 



where v\ = v^ s + v^ s ,. In this weakly coupled limit, 

S = H wc = ffiA, (6) 

in which A is the minimum of A c i = m ss / XjjVj,/\q s q s '\ or 
A qm = 2m ssl \ D v T /fi- 

Previous authors have improved the conventional the- 
ory. Lenard and BalescrP accounted for the dielectric 
nature of a plasma with a method that self-consistently 



and still diverges. Aono later showed that the Lenard- 
Balescu and Landau approaches can be combined to re- 
solve the logarithmic divergence issue, and others have 
obtained similar results with various methods F 5 " 19 ' These 
theories also obtain order unity corrections that give 
S = ln(CA), where C « 0.765. Although these self- 
consistently derive In A, they are limited to weak coupling 
because they rely on a small angle scattering expansion. 
Li and PetrasscP^ attempted to describe moderately cou- 
pled plasmas by extending the expansion procedure to 
third order; again obtaining an order unity correction. 
However, although the first two terms of this expansion 
are the only ones with C(ln A) contributions, every other 
term in the expansion has an order unity contribution. 
One must both abandon the small angle scattering ex- 
pansion, and obtain a convergent collision integral, to 
include correlation effects. 

The symmetry properties of the scattering process that 
we use are the same as those used to derive the Boltz- 
mann collision operator itself: (i) the total differential 
velocity-space volume element is invariant, d 3 vd 3 v' = 
d 3 vd 3 v' , (ii) the magnitude of the relative velocity vec- 
tor is invariant, u = u, and (iii) the "inverse" collision 



satisfies <t(v,v' — > v, v') = <t(v, v' — > v, v" 
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Here, 



<r(v,v' — > v, \r')d 3 vd 3 v' , is the scattering probability, 
which is related to the differential scattering cross sec- 
tion, er, by erdfi = / d 3 vd 3 v' <r(v, v' — > v, v'). 

Transport coefficients can be determined from velocity 
moments of the form 



(xY 



d 3 v X s(v)C B (f s J s ,). 



(7) 



Here we will be interested in Xs = ( m s, to s v, m s v 2 ), 
which are related to density, momentum, and 
energy conservation. In terms of a, Eq. Q 
is (x) s ~ s = J d 3 vd 3 v'd 3 vd 3 v' <r(v,v' -» 
t,v')ux s (v){f s (v)J s ,(v') - / s (v)jy(v')]. Making 
the interchange (v,v') «-» (v,v'), and applying prop- 
erties (i)-(iii) above, J d 3 vd 3 v' d 3 vd 3 v' <r(v,v' — » 
v, v') ux s (v)/ s (v)/ s /(v') = / d 3 vd 3 v'd 3 vd 3 v' <r(v,v' 
v,v')itx s (v)/ s (v)/ s '(v'). Thus, Eq. (7| can be written 
in terms of the distribution functions before scattering 

(X) s - S ' = J d 3 vd 3 v'{A Xs }fs(v)fAv') (8) 

in which {Ax s } = J dttauAx s , and A% s = x s (v) - 
X.(v). 

The friction force density is simply R s s — (m s v) s ~ s . 
Conservation of momentum, m s v+m s 'v' = to s v+to s /v', 
implies m s Av = to ss 'Au, so 



d 3 u {Au} I dV/ s (u + v')/ s '(v'). (9) 



The energy exchange density is Q s s = {\rn s v 2 ) — \ 
R s s . Conservation of energy, m s v 2 + m s rv' 2 — m s v 2 
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m s rv' 2 , implies Au • Au 
m ss '(v' + m ss >u/m s >) ■ Au, so 



-2u • Au, and m^Av 2 



d 3 it {Au} • I u , 



(10) 



where I u = / dV(v' - V s +m ss 'u/m s )/ s (u+v')/ s '(v')- 
Assuming the force between scattering particles is 
central and conservative, a — cr(|u|,0), so {Au} = 
-4ttmu J* d9 sin9 sin 2 (0/2) a(u, 9). Taking f„ and jy to 
be flowing Maxwellians, and assuming |V S — V a /|/t?T <C 
1, Eqs. ([9]) and ( flO 1 can then be written in the identi- 
cal form of Eqs. ^3 1 and Q, but with the generalized 
Coulomb logarithm 



2 Jo °"o 



(11) 



replacing Eq. Q for the collision frequency expression in 
Eq. (5). In Eq. ([IT]), 



4tt / d0 sin sin 2 (6>/2)cr 
Jo 



(12) 



is the momentum-transfer cross section, <j — nX^/A^ 
is a reference value for the differential scattering cross 
section, and £ = u/vt- 

If the fluid transport timescale is long compared to the 
inverse plasma frequency (^ s _ s </w ps <C 1), the polariza- 
tion response of species s is adiabatic and typically sat- 
isfies the Boltzmann relation n s = n exp(~q s (f)/T s ). As 
long as q e (j>/T s <C 1, Poisson's equation V 2 (/> = — 47r j o g 
shows that the potential around a test charge [p q = 
<?t<5(x - x t ) + ppoi] in a plasma p po \ = J2s Qs n s nas tne 
Debye screened (Yukawa) form 



-e~ r/AD , 
r 



(13) 



■x t | and XD = (Es^ 2 s r 1/2 - 



where r — 

The Debye screened potential thus requires both adi- 
abaticity (v s — s ' /u ps <C 1) and a weak interaction po- 
tential (q s (f>/T s <C 1). These conditions are often satis- 
fied in both weakly and strongly coupled plasmas. For 
electron-ion collisions, Eq. ^ shows v e ^i/ui pe ~ S/A. 
Here, A ~ nA^ has been used. For weak coupling 
(A > 1), v e ^i/ui pe ~ In A/A < 1. We will find that 
for strong coupling (A <C 1), v e -i/La pe ~ A In 2 A <C 1 [see 
Eq. (21)]. Thus, adiabaticity can be satisfied in both 
limits. For Eq. (lj), \q</>/T\ ~ (6 min /r) exp(-r/A D ). 
At the mean separation distance r ~ rT 1 ^, \qcf>/T\ ~ 
A -2 / 3 cxp(— A -1 / 3 ), which is also small for both weak 
and strong correlation. The ratio of mean interaction 
distance (n -1 / 3 ) to collision length (A e _i ~ v e -i/vTe) is 
also small in both regimes: r/A e _j ~ 5/ A 4 / 3 <C 1. 

Quantum regime. -An analytic solution of Eq. ( |11[ ) can 
be obtained in the quantum regime (A qm < A c i). Us- 
ing the Yukawa potential from Eq. ( 13 ) in the first Born 
approximation yields 



4^£ 4 [Aq„ 2 £-2 



■sin 2 (6»/2)] 5 



(14) 
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FIG. 1. Generalized Coulomb logarithm, H, from the theo- 
retical prediction of Eq. (201 (solid line), and the best fit line 
to the MD simulations of Dimonte^ (diamonds). Also shown 
are the weakly-coupled results from conventional (dotted line) 
and convergent (dashed line) kinetic theories. 



The momentum scattering cross section for Eq. ( 14 1 is 

_am 2<7 



C4 



ln(A 2 £ 2 + 1) - 



qm 



i 



A 2 £ 2 



(15) 



The effective Coulomb logarithm, 3, can be computed by 
putting Eq. (151 into Eq. (11), which yields 



1 



e A ^ J B 1 (A q ^)(l + A qi 2 ) 



1 

2' 



(16) 



where E\(x) = J°° dte t /t is the exponential integral. 
For A qm > 1, 3 qm = InAqm- (7+l)/2 + 0(A-2), where 
7 ~ 0.577 is Euler's constant. 

Classical regime.-Fov classical plasmas (A c i < A qm ), it 
is convenient to write Eq. (Ill in terms of the impact 
parameter b. Using dOsmOa = bdb, and 9 = it — 20, 
where 



9 = b 



dr 



r 2 ^\-U cS {r, b) 



(17) 



and r (b) solves U c s = 1, the momentum-transfer cross 
section is 



47T 



db b cos 8 . 



For the Yukawa potential of Eq. ( 13 ), 

A 



U, 



eff 



±- 



"A cl £ 2 r 



D e -r/A D + 



(18) 



(19) 



in which + refers to the repulsive (q s q s ' > 0), and — to 
the attractive (q s q s ' < 0) Yukawa potentials. 



An exact analytic solution of Eqs. (117|-(19) for a. 



not known, but this system has been studied numerically 
and in asymptotic limits One significant difference 
that arises for classical plasmas when using the screened 
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FIG. 2. Numerical solution of 3 c i(A c i) for the attractive and 
repulsive Yukawa potentials, along with the asymptotic so- 
lutions from Eqs. (20 1 and (211 and the quantum solution, 

1 — 'qm 

(A qm ), from Eq7"p6b. 



Coulomb, rather than bare Coulomb, potential is that 
S changes for attractive or repulsive collisions. This is 
not a feature of the quantum regime in the first Born 
approximation because the solutions merge at high en- 
ergy, where the Born approximation is valid. Hahn et 
ali^ have studied the smooth merging of quantum and 
classical solutions, which occurs when A c i£ 2 ~ A qm £, by 
accounting for the statistics of distinguishable and indis- 
tinguishable particles. 

Khrapak et alW^ have considered a s for attractive col- 
lisions between dust particles. One complication is that 
a potential barrier forms when A c i£ 2 < 1/13. Above this 
critical value, a good approximation of the momentum 
scattering cross section is <rf « 4er £- 4 ln(l + A c i^ 2 )^ 
Using this approximate cross section, Eq. (11) gives 



"Cl.c 



exp(A- 1 1 )i? 1 (A c - 1 1 ) 



(20) 



for the generalized Coulomb logarithm when A c i > 1/13, 
which covers the weak and moderate correlation regimes. 
In the weakly coupled limit A c i ^> 1, this returns the con- 
ventional Coulomb logarithm 3 c i iQ = In A c i — 7 + 0(A~ 1 1 ). 
Figure [l] shows a favorable comparison of Eq. (20) to the 
molecular dynamics simulation results of Dimontc et ali^ 
In the simulations, S was determined from thermal re- 
laxation in an electron-ion plasma: dT e /dt 2Q e ~ l /3n e . 
The green line, and triangles, in Fig. [T] represent the best 
fit line, S = ln(l + 0.7A c i), of the simulation dataP^Also 
shown are the weakly coupled results of conventional^^ 
and convergent JEHE1 kinetic theories. 

Beyond the barrier potential, A c i£ 2 < 1/13, Khrapak 
et al. show that a s ss Ax A 2 ;[ln 2 (A cl £ 2 ) - 21n(A cl £ 2 ) + 



0(1)], where A 
Eq. |TT|) gives 



0.81. In this strongly coupled limit, 



(^l/2)A 2 1 [ln 2 (A cl ) + ln(A cl )(l - 7 ) + 0(1)], (21) 



for the generalized Coulomb logarithm. Note that A/2 ~ 
2/5. The asymptotic solutions from Eq. (20) and (21) 



are shown in Fig. [2] al ong with direct numerical solutions 
of Eqs. (11) and ( 1 7 ) - ( 1 9 1 . Both the repulsive and at- 
tractive solutions asymptote to Eqs. (|20|) and (21 ) in the 



large and small A c i limits, but separate in an intermedi- 
ate region (10~ 2 < A c i < 10 2 ). Figure [2] also shows the 
approximate quantum solution from Eq. ( |16[ ). Although 
a given collision in any plasma may be quantum or clas- 
sical, depending on £, here we assume all collisions are 
either quantum or classical. Accounting for both simul- 
taneously requires 5 = S(A c i,A qm ), which may be im- 
portant in ICF plasmas that can traverse both regimes 
in a single implosion. 

In summary, the conventional calculation of fluid trans- 
port coefficients has been generalized to include large an- 
gle collisions, which are important in correlated plasmas. 
Transport is more complicated in correlated regimes, e.g., 
attractive and repulsive collisions can be distinguished, 
but the primary extension of conventional theory is a gen- 
eralized Coulomb logarithm. This theory made use of the 
Yukawa potential for individual particles, which can of- 
ten be justified regardless of the correlation strength. It 
is limited by the binary collision approximation, which 
breaks down at sufficiently strong coupling. 
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